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Abstract 

K— j We introduce a statistical agent based model to describe the phenomenon of drug abuse and its dynamical 

evolution at the individual and global level. The agents are heterogeneous with respect to their intrinsic inclination 

to drugs, to their budget attitude and social environment. The various levels of drug use were inspired by the 

professional description of the phenomenon and this permits a direct comparison with all available data. We show 

that certain elements have a great importance to start the use of drugs, for example the rare events in the personal 

experiences which permit to overcame the barrier of drug use occasionally. The analysis of how the system reacts 

to perturbations is very important to understand its key elements and it provides strategies for effective policy 

making. The present model represents the first step of a realistic description of this phenomenon and can be easily 

O generalized in various directions. 

CO 

CO 

Introduction 

CO 

The ambition to develop a quantitative description of people's behavior and introduce novel ideas and methods in socio- 
economic disciplines 1,2, 3] (see also the future ict project www.futurict.eu) is one of the main challenges of statistical 
physics and complexity theory. Agent based models (ABMs) [H [5] represent a broad framework to address these 
questions. They can describe some of the most important properties that get inspiration from physical phenomena. 
These are for example the importance of heterogeneity, large (critical) fluctuations, self-organized criticality and the 
lack of cause-effect relation. The the key features of ABM is the suitable choice of the nature of the agents, their 
interactions and their dynamical evolution. 

In this paper, we present an Agent Based Model which aims to describe, at a reasonably microscopic level, the 
phenomenon of drug abuse, its evolution and control. The scheme of the model was inspired by the participation of 
one of the authors (L.P.) to the national initiative PREVO.LAB [7J [H] whose objective is to analyze and control 
the phenomenon of drug abuse in Italy. Therefore, the model is closely aligned with the professional analysis in this 
field and permits a direct comparison of concepts and parameters with those actually observed and analyzed in reality. 
' . I In this respect, the model is rather realistic and suitable for direct applications, including control, forecast of the 
k> phenomenon and for policymaking. The information in this field is very scattered and ranges from highly accurate 
k>( information about people who become hospitalized to the little known detail about those at the beginning of the 
j_j process. The model aim is to provide a complete framework to describe the phenomenon, whose parameters are fixed 
by the best known facts. In this way, the model is able to extrapolate the knowledge of the less known submerged 
elements. This Agent Based Model, with its parameters fixed by real observations, exhibits the following features: 

• The tail fluctuations of personal experience (environment and interactions) are critical components for overcoming 
the natural barrier to start drug consumption, for the first time. This also implies that the heterogeneous nature 
of the social-network is very important and cannot be represented by an average situation. 

• The economic barrier plays a relatively minor role with an appreciable importance only at the beginning. 

• The model permits to track the individual history of an agent and in turn, this information could be directly 
compared with medical and other data. 

* Corresponding author. E-mail address : riccardo.diclemente@imtlucca.it 
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• One can analyze the response of the system to external changes of some parameters. This can help optimizing 
the control and defining a suitable policy. 

• The model is flexible and can be easily improved by introducing specific elements that may be inspired by new 
observations. For example, it could be implemented in a specific social complex network [9] and one can also 
introduce different drugs with different market socioeconomic characteristics, related for example to the age, the 
economic background or the gender. 

In summary, we introduce an ABM with variables and parameters defined in a way that makes efficient use of all 
data and information available in the field. This permits us to achieve a complete description of the phenomenon both 
at the individual and global level, and to identify its crucial elements and its responsiveness to changes in any of its 
parameters. 

The study of drug consumption and trend predictions began in the '70s with the analysis of trend historical data 
of New York City [10]. In the '90s Everingham and Rydell [11] proposed a Markovian process to analyze and predict 
cocaine consumption in the USA. The first ABM model introduced was Drug-mart by Agar and Wilson [T2][T3], revised 
by the UK Office of Science and Technology report [14) . which is a simple network model. Another more complex 
model is the Australian SimDrug |15) which simulates the consumption of heroin in Melbourne with accurate people's 
dynamic reconstruction. 

With respect to these preliminary models, our ABM has a much higher level of complexity and realism and it 
provides information on all elements of the phenomenon both at the individual and global level. To this end, it is 
important to note that our main conclusions could not have been derived from these previous models. 

Results 
The Model 

In our model, an agent i, with i G [1; N] , is characterized by a number of parameters and interacts with the environment 
and with the other agents at each time t € [1; T]. In the dynamic the unit time step will correspond to one week. We 
now briefly introduce the key elements of the model, which are described in more detail in the methods section. 

Age yi(t) — [10 : 100] years. The starting distribution of the agents' age at to is taken from the Italian data of 
ISTAT. In a first stage, we consider the property of a static situation j/j(t) = yi(t ). Later, we also introduce the 
dynamic evolution of the system by increasing the agents' age with time. 

Consumers type Si(t) = 0;1;2;3;4 defines the consumption stage and the level of dependence reached by a 
consumer: non users (Si — 0); mini user with age yi(t) < 26 (Si — 1); occasional user (Si — 2); frequent user not 
pathological (Si = 3); heavy user pathological (Si = 4). These categories are inspired by the professional classification 
of drug abuse from "Diagnostic and Statistical Manual of Mental Disorders" [16]. The mini user with yi < 26 are 
agents with low income that can be in contact with the normal dose and with a lower quantity of drug, the mini 
dose. This mini dose, as we will see later, was introduced by drug cartel with the purpose of stimulating in the young 
consumers a possible future dependency on drug [5] [7J . 

Personal budget and saving behavior. Behavioral analysis [7J has shown that the money a person would spend 
on initial drug consumption is rather related to personal saving attitude for leisure opportunities, than to the total 
personal budget. Therefore, the only variation we make in the budget rrii(yi(t)) refers to the difference between the 
adult working population with age yi(t) > 26 with m^y^t)) = 1, and young unemployed or student population with 
age yi(t) < 26 with rrii(yi(t)) = 0.2. The threshold of age 26 represents the average age of the change of spending and 
work habits [7J. The crucial point will be the parameter £ [0; oo] (saving behavior), which describes the tendency 
of an agent to save or spend money. Inspired by various indicators of social behavior, we assign to it a lognormal 
distribution [17] with variance T a and average T M . A value of 7, ~ means that the agent is inclined to spend his 
leisure budget while "ft > 1 corresponds to a strong willingness to save money and do not use it for buying drugs. 

Personal exposure ej describes every possible contact between the agent and the drug world, in particular the 
input the agent receives about drug consumption. In the present model we adopt a simplified model of interaction 
among agents but the generalization to a complex social network would be a natural extension in the model [H HH1 [22] . 

Inclination to drug fa controls the intrinsic (genetic, educational, cultural, etc.) agent's barrier towards drugs. 
Also in this case, inspired by the social studies |19j . we adopt lognormal distribution with a = B a .ji = B^. For 
every agent we define the barrier to be overcome in order to start using drugs: b*(t) = fa/a(yi(t)); where a(yi(t)) is 
a function dependent on the agent's age and it has the characteristic to increase the barrier for older agents. This 
function is described by a simple analytical form which reproduces the data of Prevo.Lab [H] . 
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Drug addiction: describes how an agent becomes addicted to drugs and how this modifies her behavior. This 
property is defined by the function di(wi(t), Si(t)) which can increase or decrease b*(t), depending on the persistence 
of the agent in a given stage, where 7Tj(£) is an exponential function of consecutive drug assumption. The inclusion of 
this effect redefines the drug dependence function: 6j(t) = b*(t) — dj(7Tj(i), Si(t)). 

In the drug market there are two different dose levels: mini-dose and dose. The first has a lower quantity and a 
cheaper price. It is used for attracting young people to drug consumption while the second represents the normal dose 
for adults and for young usual users. The mini-dose is very important for the analysis of drug consumption because 
it is often the way in which the young generations start the drug use. Moreover it represents the cheapest and easiest 
way of transmission of drug among young people [7]. In our model the mini-dose could symbolize a light drug that 
can lead to the intake of the normal one. The mini-dose price p m (t) is supposed to be substantially smaller than a 
dose. We have chosen an lower magnitude order as reasonable estimation of the mini dose price which is about one 
tenth of the dose [6]. The price of the dose Pd(t) is based on a simplified relation between supply and demand. In our 
case we assume that the total supply is fixed so the price is simply due to the total number of users. 



PS) = N Pm(t) = U 



This relation could be easily made more realistic with a variable level of total drug amount. 

The agents interact and evolve through the comparison of the two major barriers: the Economic Barrier "EB" . 
and the Socio- Emotional Barrier "SEB" . The first barrier compares the price of the drug with the buying capacity 
and it is defined by: 

f^f))<l-7, if W(*)<26 
l^OT))< 1 -^ if »(t)>26 

The second barrier has the role of describing the social opportunities and the emotional process involved in drug 
consumption. The social opportunities are defined in terms of the events by which the agents are subjected during a 
time unit. We describe this "daily noise" by ri(t), random number drawn each time for each agent from a Gaussian 
distribution, = , R a . Our choice to use the Gaussian distribution is due in large part to the lack of quantitative 
information that is available. It was also chosen because of the symmetric shape that could lead to a daily event that 
could generate a positive or negative effect in the drug intake. When more information becames available, it will be 
easy to modify accordingly. The condition to overcame the SEB is defined by: 6j(£) + n(t) < ei{t). 

We have now all the elements that define the dynamics and the interaction between the agents. In summary the 
most important parameters of the model are 7i, rj(i) and ft which are respectively money saving behavior, daily noise 
and intrinsic inclination to drug. 

We can define the transition probability of each agent from the stage <§,(<) to Si(t + I) as indicated in Figure la. 
The corresponding rates are described in detail in the methods section. 



Tuning the parameters 

The first point we consider is to tune the model parameters to a realistic description of the observed data. A crucial 
parameter in this respect is the observed fraction of people using drugs. In order to address this question we have 
made several simulations with N = 1000 and T = 1000 (the time unit t is one week). The strategy is to fix two 
distributions and examine the change of drug using population during the variation of the third one. Then we repeat 
the operation for all three distributions. In Figure lb.c.d we see the change depending on the three distributions. We 
choose to have in our population about 10% of drug users (a realistic value considering all different drugs |6, 201 l^lj. 
which could be adjusted when better data is available). When studying the real drug consumption, in addition to 
these total values, it is important to split the agents in four age ranges [5] labeled as: young (15 — 25), young-adult 
(26 — 35), adult (36 — 45) and senior (45 — 100). Also, the individual percentage of users within the four age ranges 
is used to fix our parameters. These procedures define the three distributions, which we name "Usual Distributions", 
UD: 

f ft — » tognarmal(£ / , = 1.2,B <T = 2.7) 

I 7l — » lognormal(r M =0.6,r CT = 0.9) (3) 
[ r.i(t) — » gaussian(i? Al = 0, R a = 0.6) 

Therefore, these optimized distributions provide a complete and realistic parametrization of our model which can 
now be used to investigate a variety of problems. We start by considering the individual personal history per agent. 
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Fi gUre 1'. Stage transitions and parameter settings, a. All possible stage transitions for an agent at stage Si(t). b. c. d. Simulations 
intended to select the main parameter of the model in relation to the realistic stage of the phenomenon (see text for more details). In particular 
in b. we show the effect of the parameter /3i (inclination to drug) distribution (S m ;B CT ), in c. effect of a parameter ji (saving behavior) 
distribution (r m ; r CT ) and in d. Ti (daily noise) distribution (K m — 0; R&). The blue circles identify the realistic areas to fix the values of the three 
distributions so that they correspond to a total of 10% of users. This parameter optimization refers also to the age ranges discussed in the text. 



Individual agent history 

It is interesting to note that the individual history can be vastly different in reality and we are going to observe that 
our ABM is able to reproduce a dynamic variety of agent's history. The table in Fig. 2a shows At = 68 consecutive 
iterations of eight users, and it outlines some of the different users typology that can be produced by the model: 

• #64 is a senior agent that never tested drugs; 

• #141 is a young agent that, as an occasional consumer, tries the normal dose; 

• #l(young-adult) and #343(young) are two agents that are occasional consumers; 

• #15(young-adult) and #130(young) are both heavy consumers that evolve into addicted; 

• #45(adult) and #494(young) are two agents that try drugs on extremely rare occasions. 

The fact that our model reproduces a realistic heterogeneity of agents' history is a very important element. This 
realistic variability is in fact crucial to have an accurate description of the phenomenon that could be compared to a 
considerable amount of real data. A possible upgrade would be analyzing and contrasting the real drug user clinical 
histories. By means of the above we could provide a probabilistic prediction of detoxification options for each clinical 
case. In order to gain more insight in the individual evolution, we can observe the trend of the barrier for inclination 
to drugs use, the Social Emotional Barrier (SEB) and the stage reached by an individual user i at anytime t. The 
Figure 2b shows the history of a young user yi = 19 with a low intrinsic inclination barrier fii, for At — 300 consecutive 
iterations. The user stage Si(t) is indicated in red dots, the left term of SEB inequality barrier in green, and in blue 
the starting value of inclination to drug b*(t). We can see that this agent develops all the stages with important 
fluctuations and finally ends up being an addict. At this level the green barrier line declines strongly and finally we 
can also observe that the agent detoxifies. In the case above the effect of aging is not yet included b*(t) = b*(t ). 
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Now we consider a complete dynamic evolution in which aging is explicitly included. We set the year length 
r = 50t, in which one iteration corresponds to a real week time spend. We elaborate a set of rules to define the death 
of an agent due to the aging and the birth of the new one in order to have a constant number of agents N. We also 
define a simple process, from the data [51 [TSl HO]: which describes the possibility that an agent can die of overdose 
and be replaced by a new one. In Fig. 2c we reproduce the personal history of an agent in the aging process. The 
figure shows the increase of the barrier to b*(t) due to the agent's aging. The agent starts with a high value of /3, 
(low tendency to use drug) but the daily noise induces to try several drugs (Si(t) = 2). However, when the agent 
grows up, b* (t) rises, so the daily noise becomes less effective and finally the agent no longer uses drugs. This example 
shows that the possibility to monitor a single agent's history can be an important tool to compare with real medical 
data. Once again the future analysis could study the probability of stage change and compare it with real data from 
detoxification centers. 
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Fi glire 2; Individual agent histories, a. The tabic shows the histories of eight agents for At — 68 consecutive times (One time unit corresponds 
to one week). One can see that our model is able to reproduce a variety of a realistic heterogeneity in agents' history. This history ranges: from the 
agent #64, who never tried drugs, to agents characterized as heavy users #15; #130. b. Personal agent history without aging: the agent becomes 
a heavy user because after various fluctuations she reaches the stage Si(t) — 4. The green line describes the left term of Social Emotional Barrier 
and one can see that this barrier is drastically lower when the stage of dependence is reached, c. Personal agent history with aging. In this case 
the agent is a casual user, meaning that agent has frequently used drugs at a younger age. Then the barrier of inclination to drug use (blue line) 
b*(yi(t)) increases with the age. We can see that, in this way, the use of drugs Si(t) — 2 becomes more and more rare. Finally when the agent gets 
older the use of drug is completely eliminated. This example shows that the model provides a detailed description of individual histories, which 
could be directly compared to real data. 



Global property and effects of perturbations 

Our model also allows us to examine how the users population reacts to various sources stress, ranging from a different 
opportunity to be involved in drug use (r^), a new agent's intrinsic education (fli) or a change in saving behavior 
(ji)- With the purpose of analyzing the agent's dynamic under these social modifications, we fixed two of the UD and 
analyzed the effect of changing in-time of the third one. We will activate the in-time modification between the time 
steps, t + < t < 2t + . 

In Figure 3a we change the daily noise distribution from the UD to = 0.3, R a — 1.6, with t + = 4r. The Figure 
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shows how an increase of noise leads to an increase in the number of agents in the drug stage. As it is shown from the 
insert of Figure 3a the ri(t) distribution has different impact for each age range, given the dependency of b*(t) on age. 

Figure 3b exhibits the change of the population stages Si(t) due to the modification of the intrinsic inclination 
toward drug use distribution (B^ = 1.2, B„ = 2.7) for the new born agent between 20r < t < 40r. The agents who are 
born during the perturbation periods are more inclined to use drugs than the others. The insert of Figure 3b shows 
the trends between the age ranges and how the aging of the new agents increases temporally the number of drug users. 

In Figure 3c we change the tendency of agent spending behavior. We replace the distribution of 7, from the UD 
with (r m = 0.1; T a — 0.3) for the new born between 20r < t < 40r. During this perturbation period the new born 
agents are more inclined to spend money for leisure, in this case drug. One can see from the insert of the picture 
the different trend among age ranges. According to the analysis, the saving behavior does not appear to have an 
important role in explaining drug assumption. The decrease in-time the % distribution produces a limited impact in 
the number of users compared to the effect /3j. 

Finally we focus on the role that mini users have on the introduction of drug among the new generations. This 
is an important point because it is hard to monitor and as an appreciable effect only after a long time. In Figure 
4, we compare two situations of a population evolution. For the first 15r years the simulation does not consider the 
mini-dose. Consecutively we included the mini user. We can see that the existence of mini user stage leads after many 
years to a large increase in the total number of drug users. 



Discussion 



In conclusion, we have developed an agent based model for the phenomena of drug use, which is relatively simple, 
but at the same time rather realistic. Its concept and parameters have been established in close relation with the 
professional studies of the real phenomenon. The model permits the study of the problem at the individual and global 
level. The individual histories of the agents can be extremely diverse and appear to provide a realistic description of 
the real situation which could be easily compared to medical data directly. At the global level we can identify the 
crucial most important parameters, as well as those that play a minor role. This leads to a basic understanding of 
the complexity of the phenomenon and to the possibility of analyzing how the situation changes if the parameters are 
modified. This point is essential in defining the strategies for controlling and reducing the phenomenon. 

One of the main conclusions of our studies is the fact that drug use is usually trigged by a rare event in the 
personal experience. This shows the importance of the tail of these heterogeneous experiences, which allows us to 
predict that the introduction of additional heterogeneity in the structure of the society will lead to important effects. 
In particular, our model's interactions among agents are assumed to be similar on average. Therefore we expect that 
the introduction of a complex network distribution of this iteration will be relevant. In this perspective we can already 
predict that the hubs of this distribution will be very important for policymaking. 

The model is meant to represent a realistic basis of analysis and discussion, and can be easily generalized with a 
complex network of social interactions and with multiple drugs with different social characteristics. We believe that 
these studies can lead to a higher level of understanding of these phenomenona and can be useful for a more effective 
policy making. 



Methods 

Model parameters specification 

Following we explain in detail the parameters that characterize the agents. 



Personal exposure 

The personal exposure e^, describes every possible contact between the agent and the drug world and in particular 
describes how agents feel about drugs in general and their consumption in particular. It is defined by: 



i(t) = 



N 

J2 6(S k (t)=l, 2,3,4) - 
k=l . 

N -1 



5(S i (t) = l,2,3,4) 



(4) 



This expression simply counts the fraction of drug users within the entire population, but does not include the single 
agent under observation. The latter represents a type of mean field interaction, which could be easily generalized to 
complex network situations [20 [22] . 
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Fi glire 3'. Response of the system to external changes, a. Fluctuation of the life style. Here we test how the system reacts to increased 
fluctuations for the rare events (tail of the ri(t) distribution) that induce the starting of abuse. During the period 4r < t < 8r the distribution 
of daily noise increased (R m — 0.1; -R CT — 1.3), the various ranges of population arc affected by this change. In the insert we show how the total 
number of drug users is distributed among the various age ranges. From these studies, one can conclude that the system is strongly sensitive to this 
tail fluctuations, which represents a very important point for the control of this phenomenon. b. Cultural fluctuations. Here we show the effect 
of the changes in the genetic cultural barrier of inclination to drugs. Wc change the distribution of inclination to drug 04 for the new agent during 
the period 20r < t < 40r with (B m — 1.2; B a — 2.7). The lower intrinsic barrier leads to a creation of a new users generation more inclined to 
drugs abuse. We can see that the effect is not instantaneous and wc can follow how it develops over the years. Also in this case, as we show in the 
insert, one can see the large difference of age effect between the various age ranges. The pattern revealed can be directly compared to real data and 
is important for monitoring and controlling the phenomenon, c. Fluctuations of saving money behavior. Here we examine the effect of the change 
of the saving behavior. Wc change the distribution of -yi for the new agent during the period 20r < t < 40r with (r m — 0.1; T CT — 0.3). A lower 
value of 7^ means that an agent tends to spend money easily. In particular a value near to zero for an agent means that she has no problems for 
drug purchase. Wc can sec that the effect is similar to the cultural fluctuation but it is of minor impact. The effect of spending capacity decreases 
with the age and in particular it does not affect the senior agents. 
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Figure 4 '. Mini dose activation. Wc consider a population of N — 7000 agents. At the beginning the market is without the mini dose, at a 
time t — 15r the mini dose is introduced in the system for young people. We only show the drug user stages and one can see an immediate change 
of the global users followed by a further increase of addict population. The results show that the effect of mini dose for young generation represents 
and extremely important clement for the dynamics of the system. 



Age function 

The function a(yi(t)) is dependent on the agent's age and its purpose is to increase the barrier for older agents: 

a{ Vi {t)) = ^ L (5) 

The values of a = 1; d = 0.5; c = 3 are chosen in order to outline the importance of age in the consumption of drugs 
[7] . This expression is a simple mathematical representation of how age influences the probability of drug assumption 
Fig. 5a. The main characteristic of this function is that the peack of this function is located between sixteen and 
twenty-four years and then it slowly decreases, as it has been shown by the studies of Prevo.Lab. As we can see from 
the Figure 5b the a(yi(t)) increases the average value of the b* due to age behavior. 

Drug addiction 

The drug addiction describes how an agent becomes addicted to drugs. This property is defined by the function 
dependence di(iTi(t), Si(t)) which has the purpose to decrease the barrier b*(t). The function 7T;(t) "stage permanence" 
is an experience counter: 

7^ = 1 if S t {t)*S t (t-l) 

7T i (t)=7r i (t-l) + l if Si(t) = Si(t - 1) W 

and the function dependence is: 

difa(t),Si(t)) = I(t)(exp(7r i (t)5 i (t)*x(2))-l) + di(7r i (t-l),S i (t-l)) (7) 

This expression represents an exponential growth of the addiction parameter with the persistence of the certain 
drug user stage. The term x( z ) represents the drug coefficent, which describes how addictive the drug of type z is. In 
this regard, it would be easy to introduce drugs of different types. For the present study we consider a single case in 
which x( z ) — -0001. This parameter could be changed for increasing or decreasing the effect of the drug addiction. 
The parameter li(t) is the first experience coefficient, which describes the quality of the agent's first experience in 
using the given drug: 

k(t) = -2 if ki(t) > 0, 75 & Si(t) < 3 . . 

k(t) = 1 Otherwise [ ' 
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Figure 5: studies of barrier age dependence, a. The function age coefficient (a(yi)) is described by a simple analytical form which is a 
reproduction of Prevo. Lab's data. b. This is an example of a distribution of the intrinsic agent barrier with age dependence b*(yi) for a simulation 
of 10000 agents. 



where k t (t) is a random number with uniform distribution between [0:1]. The two possible values of Z$ (t) are chosen 
to give a mathematical representation of the social hypothesis, that a negative experience has more psychological 
impact than an equivalent positive one |23j . We can redefine the barrier of inclination to drugs with 6j(t) = b*(t) — 
di(iti{t), Si(t)) in order to outline the role of dependence and first experience in the drug consumption and addiction. 

Rules for the stage transitions 

We define the rules that enable agents to have stage transitions from stage Si(t) to Si(t + 1). The Figure la in the 
paper shows all the possible transition stages. Each transition has some different condition due to the initial stage 
Si(t). 

From Si(t) = 

From the stage Si(t) = we can have stage transitions into: 

• Si(t + 1) = 1 if the EB and SEB barrier are verified for an agent with age Ui(t) < 26; 

• Si(t + 1) = 2 if the EB and SEB barrier are verified for an agent with age yi(t) > 26; 

• otherwise the agent remains in the stage of non user Si(t) = 0. 

As we can see from Eq|9]the difference between the final stage reached by the agent is dependent only on age: 

Si(t + l) = l if [ Vi < 26 & < 1 - 7l & bi(t) + n(t) < ei (t)} 

Si(t + l) = 2 if [ yi >26k ^fgj <l-7i &6 i (i)+r<(t)<e i (t)] (9) 
5i(t + l)=0 Otherwise 

From Si(t) = 1 

Only an agent with yt(t) < 26 can reach this stage. The possible transitions as we show in EqflO|are: 

• Becomes non-user, if EB or SEB are not verified; 

• Persists in the stage of mini-user if the agent gets the sufficient daily noise (opportunity) and has a budget; 

• Change in Si(t + 1) = 2 if the agent has funds to purchase the full dose and the SEB is verified without need 
of the daily noise. The possibility for an agent to get in contact with the full dose is dependent on personal 
behavior associated with dependence level rather than the daily occasions. 
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Si(t + 1) = 2 if \p d (t) mi {yi) < 1 - 7< & h{t) < «*(*)] 

Si(t + 1)=Q if [p m {t)m t (y t ) > 1 - 7i or + r,-(t) '> e;(i)] (10) 

£j(t + l) = l Otherwise 

An agent becomes Si(t + 1) = 2 when the function dependence and the first experience coefficient decreases the 
SEB. 

From Si(t) = 2 

In this stage the agent is an occasional user of normal dose, the agent can: 

• Return to non-user stage, if EB or SEB are not verified; 

• Stay in Si(t + 1) = 2 if both conditions are verified; 

• Become a frequent user if the agent's inclination barrier to drug is reduced by the dependence (SEB without 
daily noise). 

The following equation shows the transition possibilities: 

' S,(f + 1)=0 if [^g)>l-7, or b i (t)+r i (t)>e i (t)] 

< S,(t + 1)=3 if [^g y <l_ 7i & b i (t)<e i (t)} (11) 
Si(t + 1) = 2 Otherwise 

As we can see, the transition rules are quite similar to the previous stage. We have considered that the stage Si(t) = 2 
is reachable by all agents without being contingent on age. 

From Si(t) = 3 

The agent is a frequent user. She does not care about the EB because the addiction raises the necessity of a dose and 
the money becomes a negligible constraint. She can become: 

• Heavy user, pathological, if her inclination barrier to drug decreases by a factor 5 or more; 

• She can decide to stop using drugs if the social environment changes drastically, or with the 10% probability, 
which represents the possibility to stop using drugs due to exogenous environmental factors. 

• Otherwise the agent i continues to take drug as a frequent user. 
Here we show the transition rules: 

Si(t + 1)=4 if 5b t (t)<e z {t) 

Si(t + 1) = if [bi(t) > e t (t) or a(t) > 0.9] (12) 
Si(t+1) = 3 Otherwise 

Where Ci(t) is a random number from uniform distribution [0 : 1]. There is no more the daily noise f"i(t) and the 
transitions depend on the variation of bi(t) due to the dependence function and the stage persistence. The factor 5 of 
the first condition is based on our decision [SJ [7J [5] to discriminate in a significant way the difference between the two 
stages Si(t) — 3 and Si(t + 1) = 4. The bi(t) > ej(i) therefore represents the decreasing possibility of the exposure 
factor due to a sudden change in the social interactions. 

From Si(t) = 4 

When an agent becomes a heavy user, she has only two possible choices according to the transition table showed in 
the main article. The possibilities are: 

• The agent can leave drug with the probability of 5%; 

• The agent can persist in the stage Si(t + 1) = 4. 
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At this stage, the situation of the agent is independent of the condition EB and SEB due to every addiction. Therefore 
we assume only a simple probabilistic transition: 



S<(i + 1)=0 if k(t)>0.95 

Si(i + 1)=4 Otherwise { ' 

Where U (t) is a random number from uniform distribution [0 : 1] . The probability of the 5% is an estimation suggested 
by the Prevo.lab studies and is related to the possibility of recovering from drug addiction. 

In this way we have defined all the rules that allow an agent to evolve and reach all possible stages of consumption. 
So far, the agent cannot die and the age yi(t) is extracted at the beginning of the simulation and remains constant for 
every t. Later we define the rules that describe the death and birth of the agent and its aging dynamics. 

Aging 

Here we show the rules that define the process of an agent's vital cycle (death and birth), and the dynamic that leads 
an agent to overdose. 

Death and birth condition 

Based on the 2007 ISTAT data we consider the probability of death in a given year for a person of a given age. We 
label this probability P yt (t) and every r = 50i consecutive iteration we check this rule for each of N agents: 

U M < P vi{t) V& + !) = V&) + 1 (u) 
Otherwise i dies yi(t + 1) = 10 

Where Ui(t) is a random number from uniform distribution [0 : 1]. When the agent i dies, she will be replaced by 
another agent with age j^(f + 1) = 10 and she will be characterized by the new extraction of the parameters /3j,7i 
from the current value of the distributions at time t. Instead, if she lives, she will be one year older. 

Overdose rule 

By adding the possibility of death, it is possible to improve also some rules of the Si(t) = 4 stage transition. An agent 
in Si(t) — 4 can become a non-user by either detoxifying in a treatment center or by dying of overdose. When the 
agent dies of an overdose, the agent will be replaced by a new agent with the new extraction of the parameters /?i,7i 
from the current value of the distributions at time t. The following equation expresses the overdose death condition: 

Si(t + 1)-»0 if fft(«) > 0.1 lives , . 

Si(t + 1)->0 Otherwise dies ^ ' 

Where gi(t) is a random number from uniform distribution [0 : 1]. The agent with the probability of 10% dies. 
Otherwise she lives. The probability to survive after an overdose is taken from the data of SimDrug and Sert [HI [H] • 
With our dynamics we can replicate the same overdoses trend of SimDrug model too. 
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